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We construct a Galilean invariant low-energy effective field theory of boson-fermion mixtures 
and study bound fermion states on a vortex of boson superfluid. We derive a simple criterion to 
determine for which values of the fermion angular momentum I there exist an infinite number of 
bound energy levels. We apply our formalism to two boson-fermion mixed systems: the dilute 
solution of 3 He in 4 He superfluid and the cold polarized Fermi gas on the BEC side of the "splitting 
point." For the 3 He- 4 He mixture, we determine parameters of the effective theory from experimental 
data as functions of pressure. We predict that infinitely many bound 3 He states on a superfluid 
vortex with I = — 2, —1,0 are realized in a whole range of pressure 0-20 atm, where experimental 
data are available. As for the cold polarized Fermi gas, while only S-wave (I = 0) and P-wave 
(/ = ±1) bound fermion states are possible in the BEC limit, those with higher negative angular 
momentum become available as one moves away from the BEC limit. 

PACS numbers: 03.75.Lm, 03.75.Ss, 67.60.Fp, 03.65. Ge 



I. INTRODUCTION 



fermion dispersion curve e(p) sits at zero momentum: 
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Quantum degenerate mixtures of bosons and fermions 
provide interesting playgrounds where effects of quantum 
statistics become explicit. Several experiments on the 
boson-fermion mixtures have been performed, first on the 
classic example of the 3 He- 4 He mixture Q, Q and more 
recently in cold atomic gases such as 6 Li- 7 Li 
23 Na H, 40 K- 87 Rb 0], and metastable triplet 3 He*- 4 He* 
mixtures Q. Molecular condensates, formed through 
the BCS-BEC crossover, provide another possibility for 
the boson-fermion mixtures S S E3> El 0- Stimu- 
lated by experiments, theoretical studies of the boson- 
fermion mixtures predict new phenomena such as phase 
separation between Fermi gas and Bose-Einstein conden- 
sates 13] or sympathetic cooling of Fermi gas by the 
superfluidity of Bose-Einstein condensates [Tj| . 

The main goal of this paper is to investigate the bound 
states of fermions on a superfluid vortex. Superfluid vor- 
tices were realized by the recent experiment in cold Fermi 
gas throughout the BCS and BEC regimes 15] . The 



e(p) = e 



2m* 



(1) 



problem has been considered in the literature |16tll7llT8| . 
with most of the treatments following the original ap- 
proach by Caroli, de Gennes, and Matricon 0] (or 
the more rigorous approach based on the Bogoliubov- 
de Gennes equation |2(J). The reason for us to revisit 
the problem is that the Caroli-de Gennes-Matricon for- 
malism is expected to work well only in the BCS regime 
where there exist well-defined quasiparticles. Its applica- 
tion to the BEC regime where two fermions form a bound 
molecule would require an exceedingly large number of 
mean-field quasiparticle states for an accurate descrip- 
tion of the strongly bound molecules. It also cannot be 
applied to the important case of the 3 He- 4 He mixture. 

We are interested in the case when the minimum of the 



where m* is positive and quartic and higher order correc- 
tions are negligible. This case is opposite to the situation 
in the BCS theory, where the minimum is near the Fermi 
momentum, but it holds for the solution of 3 He in 4 He 
and also for cold atom gases sufficiently deep in the BEC 
regime. Our technique relies on the use of an effective 
field theory. The effective field theory formalism is es- 
pecially powerful in strongly coupled systems where per- 
turbative methods are not applicable. In our cases the 
details of the interactions become unimportant, as they 
are encoded into a few low-energy parameters. Thus the 
theory developed in this paper is universal and can be 
applied to many different physical systems. 

In writing down the effective field theory, we put spe- 
cial emphasis on the Galilean invariance [2l], |22, |23j . The 
Galilean invariance as well as the global symmetries of 
the system considerably restrict the possible form of the 
effective Lagrangian. We show that coupling constants 
appearing in the effective field theory can be expressed 
through the effective fermion mass and the derivative of 
fermion energy gap with the chemical potential. The lat- 
ter is directly related to the Bardeen-Baym-Pines (BBP) 
parameter. Both the effective fermion mass and the BBP 
parameter have been measured by experiments for the 
3 He- 4 He mixture. 

After the construction of the effective field theory, we 
use it to the study of bound fermion states on a vor- 
tex of the boson superfluid. We derive a simple criterion 
which tells us at which values of the angular momentum 
I there are an infinite number of bound energy levels. We 
will apply our effective field theory to two boson-fermion 
mixed systems, a dilute solution of 3 He in 4 He super- 
fluid in Sec. HTland the BEC regime of the cold polarized 
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Fermi gas in Sec. lIIII in order to predict possible angular 
momenta for bound fcrmion states. 



II. DILUTE FERMIONS IN BOSON 
SUPERFLUID 

A. Galilean invariance and effective field theory 

Here we demonstrate how to construct an effective field 
theory of dilute fcrmions in a boson superfluid in accor- 
dance with the Galilean invariance and global symme- 
tries. Since the microscopic theory is Galilean invariant, 
its low-energy effective Lagrangian should satisfy the fol- 
lowing identity required by the Galilean invariance: 



T 0i = mJ. t 



(2) 



where Tgi is the momentum density, m is the mass of the 
particle, and Jj is the particle number current pH |22|. 
If there are more than one species of particles in the sys- 
tem, the right hand side of the equation is a sum over all 
species. This identity simply indicates that the momen- 
tum density should be equal to the mass carried by the 
particle number current in the Galilean invariant system. 
The Galilean invariance as well as the global symmetries 
of the system considerably restrict the possible form of 
the effective Lagrangian. For definiteness, we consider 
the particular case of the dilute solution of 3 He in the 
4 He superfluid from now on. However, the universality 
of the effective field theory means that our discussion, 
with only minimal modifications, will be applicable to 
any other dilute fermion excitations in a boson super- 
fluid. 



1. 4 He superfluid phonon 

First of all, let us consider an effective field theory for 
the pure 4 He superfluid. The physical degree of freedom 
at long distance is the phase of the condensate d, which is 
defined as (*) = |(*)|e w with f being a field for the 4 He 
atoms. Under the U(l) symmetry associated with the 
number conservation of 4 He atoms, the field $ transforms 
as i) — > i9 + x- Therefore the effective Lagrangian of 
the 4 He superfluid obeying this symmetry should be a 
function of derivatives of the field # and did. We note 
that terms where derivatives act more than once on one 
field such as d or d 2, d are also possible. However, they 
are in higher orders in the power counting scheme where 
d and did are regarded as order O(l). In consideration 
of the rotational symmetry, the leading order effective 
Lagrangian £4 should be written by a polynomial of d 
and (dd) 2 as follows: 



£ 4 = £ 4 (tf, (dd) 2 ). 



(4) 



Now we impose the Galilean invariance Toi — 7714 J 
on the effective Lagrangian, where 7714 and jf 5 are the 



mass and the number current of the 4 He atoms. Since 
the momentum density and the particle number current 
are given by 



dd 

j( 4 ) = -^1 = 2W- 5U 



8d,d 



s(dd) 2 ' 

(4) 



the Galilean invariance T$i = m^Jl results in 
8C4, _ 8C4 

sd ~ 21714 smy 



(4) 
(5) 

(G) 



This equality requires that the effective Lagrangian de- 
pends on the field d only through the combination 
d + (dd) 2 /2m 4 as follows: 



£4 = P -0 - 



(dd) 2 

2m4 



(7) 



Here P(- ■ ■) is an arbitrary polynomial, which will be 
identified to the pressure of the pure 4 He superfluid as a 
function of the chemical potential later. 

In the superfluid phase, the symmetry d — > d + x is 
spontaneously broken. At finite chemical potential /J4, 
the ground state of the superfluid system corresponds 
to do — —fi4t 22]. Then we expand the field around 
the ground state as d = do + (p, where fluctuations of 
ip around zero corresponds to superfluid phonon excita- 
tions. Substitution of the expression i9 = — fi^t + if into 
Eq. results in 



£4 = P /i 4 - (p 



{dipf 



(8) 



Now we can show that the function P( - ■ ■ ) is identical to 
the pressure as a function of ^4 at zero temperature up 
to an irrelevant constant. For that purpose, we calculate 
the number density of 4 He atoms 114 by differentiate the 
Lagrangian with /14 and we obtain 



n 4 (/U4 



<9£ 4 
9^4 



dP 

(9^4 



(9) 



This equation implies that P is a function of the chem- 
ical potential satisfying P'(^i) = ^4(^4), which means 
that P(ijla) is identical to the thermodynamic pressure 
up to an irrelevant constant. Once equation of state of 
the 4 He superfluid -P(/i4) is given as a function of the 
chemical potential, the low-energy effective field theory 
of superfluid phonons is simply given by replacing [14 with 

/I4 — If — (d(fi) 2 / 2-1714,. 

In order to proceed our analysis further, we expand 
Eq. JSJ) up to the second order in fields 



(3) £ 4 =a P(/m) - 



3P 
dfi4 



{dyf 



P(H4 



«4<^ + 3 7T 



dri4 (f 
8^4 2 



2m4 
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ld 2 P 



2dfj, 2 

n4 (df) 2 

7714 2 



if- 



2m4 



(10) 
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The first term gives the pressure without phonon excita- 
tions. The second term is a total derivative of the field, 
which does not affect the equation of motion. The third 
and fourth terms represent the propagation of phonon 
with its sound velocity [22, 01 



7774 drii 



dP 



7714 dri4 



(11) 



The higher order terms which are not shown in Eq. i|10|) 
represent self-interactions among phonons. 



2. Minimal coupling between He and He 

Next, we consider the coupling between the superfluid 
phonon ip and a 3 He atom ip. One coupling term can 
be written down from the following argument. We note 
that there is an energy cost to introduce a single 3 He 
atom into the He superfluid. This energy is some func- 
tion A(/i4) of the 4 He chemical potential /J4. However, 
Galilean invariance tells us that /J4 always enters the La- 
grangian in the combination /J4 — ip — (dip) 2 j2m^. Thus, 
the Lagrangian contains the following term: 



r — 

''-gap — 



-A /i 4 - (f> - 



(9<P) 2 

2mn 



dA 

-A( fM )^iP+ — ¥>• 
C7/4 



ip'ip 
2rri4 



t^ip. (12) 



The first term represents the energy cost of introducing a 
3 He atom into the pure 4 He superfluid. The second term 
proportional to <9A/<9^ 4 gives a Galilean invariant cou- 
pling between the superfluid phonon and the 3 He atom 
put into the 4 He superfluid. However, this is not the only 
coupling between 3 He and 4 He, as we shall see below. 



3. 3 He kinetic term 

Now, let us consider the kinetic term of the 3 He field 
ip. The 3 He atom put into the 4 He superfluid has the 
effective mass 7773, not equal to the bare mass 7773 due to 
the strong interaction with the 4 He superfluid. Thus the 
kinetic term of the 3 He atom is 



restore the Galilean invariance. Therefore, the Galilean 
invariance imposes some relationships between the inter- 
action terms and the 3 He effective mass. This situation 
here is reminiscent of the Fermi liquid theory, where there 
exists a relationship between an effective fermion mass 
and a Landau parameter. 

One way to construct the Lagrangian that obeys the 
Galilean invariance is as follows [2l|. The effective mass 
of 3 He atom 7773 originates in the fact that the 3 He 
quasiparticle in the 4 He superfluid entrains superfluid 
4 He atoms due to the strong interaction between them. 
Therefore, an elementary excitation of the system, or the 
3 He quasiparticle ip, should be written as ip = e %rw ip. 
Here 77 is a parameter defined by mj = 7713 + rj 7774 , which 
represents a "fraction" of superfluid 4 He atoms in the 3 He 
quasiparticle. Using the 3 He quasiparticle field ip instead 
of the bare 3 He field ip in Eq. . we have 



iip^doip \dipY 



2mi 



rj dip 



iip* dip 
2ml 



■ HZ^Ip 



rj(p + rj' 



2ml 



(14) 



ip' ip . 



We can verify that this modified Lagrangian Ijl4(l obeys 
the Galilean invariance. The momentum density of the 
system is given by 



rr ^ 3 a / 1 <^3 q ,+ , <^3 Q 

Tot = — y diil) + —7- d t ipi + — diip 
Sip Sip^ oip 



iip^ diip 



rj diip ip^ip . 



(15) 



On the other hand, the number currents associated with 
the 3 He atoms and 4 He atoms are, respectively, given by 



SC 3 

Sdiip 



iip - 



8C3 
5di^ 



iip^diip 



2ml 



and 



sc 3 

Sdiip 



n- 



iip* diip 
2ml 



r/ 2 ^ip f ip. 



dif /t , 
777,3 

(16) 



(17) 



U - 



iip^doip \dipY 



2ml 



+ fi 3 ^ip (13) 



with fj, 3 being the chemical potential for the 3 He atom. 
However, this Lagrangian does not satisfy the Galilean 

(3) (3) 

invariance = m 3 J^ , where J, is the number cur- 
rent of the 3 He atom. This is because the Galilean invari- 
ance condition involves the bare mass, while the kinetic 
term in the Lagrangian involves the effective mass. 

The resolution to this apparent paradox is that there 
are interaction terms that contribute to both the momen- 
tum density and the particle number current in order to 



Recalling the definition of 77, the Galilean invariance 
To; = 777,3 + is indeed satisfied by the La- 

grangian l|14f> . We note that the 3 He self- interactions are 
negligible in the dilute limit of the 3 He density, while 
in general they appear as higher order terms in the La- 
grangian. 



4- Effective Lagrangian 

Finally, getting Eqs. ifTU)) . lfT2|) and (jTI)) together, we 
have the low-energy effective Lagrangian of the 4 He su- 
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perfluid phonons and 3 He excitations as follows: 



£ e g = 



2 V 2 



g x dip 



(9^r 

2m% 



2w 4 



+ (n 3 - A) ^ 



(18) 



52<^> + 53 



(dip) 



2 i 



2m 4 



ip*tp . 



Here we have introduced low-energy parameters as 



dri4 

dm 



f = n± 

m 4 



and couplings as 

1714 



5i = V- 



92 



OA 



V, 53 



OA 
<9/j 4 



>1 



, m 4 



(19) 



(20) 



While there are three independent terms representing in- 
teractions between the 4 He superfluid phonons and 3 He 
excitations, their couplings are not independent. We 
have one constraint on the three couplings 

m 3 



53 = 52 



m 4 



-5i 



(21) 



as a consequence of the Galilean invariance. Moreover g\ 
is completely determined by the 3 He effective mass m\, 
and 52 (and therefore g-$) is determined by m?. and the 
derivative of the energy cost for introducing a 3 He atom 
with the chemical potential A'(/i 4 ). These free parame- 
ters should be either determined by experiments or com- 
puted from microscopic theories. We emphasize again 
that our low-energy effective Lagrangian (|18[) and its con- 
sequences are applicable to any other dilute fermion ex- 
citations in a boson superfluid. Details of strong interac- 
tions among bare particles are encoded into the effective 
mass to* and the energy cost function A(/i). 



B. Bound states on a superfluid vortex 

One of the consequences derived from the effective La- 
grangian is on bound 3 He states on a vortex of 4 He su- 
perfluid. Let us consider a single vortex with its winding 
number w in the 4 He superfluid. The phase ip around the 
vortex is given by (p(t, r) = w9 in cylindrical coordinates. 

The equation of motion of the field ip from Eq. i|18|) 
gives a Schrodinger equation for a 3 He atom on the vortex 
ip = w6 as follows: 



Eip(r) 



d 2 ip(r) igiw dip(r) g 3 w 2 



2toj 



to 4 t 



89 



2m 4 r' 



-^(r). (22) 



E is the energy eigenvalue of the 3 He atom and is negative 
for bound states. Separation of variables leads to the 
equation for the radial wave function R(r) as 



(2mlE-kl) 



dr 2 



R(r) 

ia_ 

r dr 



251 



wl 
m 4 r 2 



53 j 

m 4 r A 



(23) 
R(r). 



Here / is an angular momentum of the 3 He atom and 
k z is its momentum along the vortex line. Introducing 
E = E- k 2 j2ml < and 



53— w 2 
rriA 



we can rewrite Eq. lll'.'il) as 

8 2 R 



ml 
2g 1 —±wl 
m 4 



dr 2 



ldR 

r dr 



2ml \E\ 



R = 0. 



(24) 



(25) 



Solutions of this equation are given in terms of the mod- 
ified Bessel functions as R(r) = I±i K (\j2m% \E\r). 

One should remind that our effective field theory writ- 
ten in terms of the phase of the condensate is valid only 
far away from the vortex core where the magnitude of the 
condensate is almost constant. Suppose the Schrodinger 
equation l|25l) is valid for r > r$ and the binding energy 
IE 1 1 is small enough to satisfy \j2ml \E\ r -C 1. In this 
instance, the solution of Eq. I|25|) near the point r = r*o 
turns out to be 



R(r) = Csin 



Kin \ 2ml \E\ 



(26) 



with C and <fi being arbitrary constants. This radial wave 
function for r > ro should smoothly connect with the 
radial wave function from r < ro. If the binding energy 
\E\ is small enough compared to the potential energy 
at r < ro, the radial wave function for r < ro will not 
depend on \E\. Therefore, the solution l|26() is required 
to be independent of \E\ at r = r . 

The logarithmic derivative of Eq. (|26|l at r = r results 



R'jro) 
R(ro) 



= — cot 

ro 



Kin \ 2mt \E\r 



(27) 



In order to satisfy the requirement, allowed energy levels 
\E\ should be discretized as 



n\n(^2ml \E n \ r ) 



D — 717T 



or, equivalently, 



E,, — 



e 2D/ K 

2m*r 2 ' 



-2nir j k 



(28) 



(29) 



Here D is an arbitrary constant and n is an arbitrary in- 
teger large enough to satisfy \E n \ <§; l/(2m3r 2 ). There- 
fore, as long as k 2 is positive, we have an infinite number 
of energy levels for bound 3 He states. While we can not 
determine the absolute values of bound energy levels in 
our approach, their ratios are independent of unknown 
constants ro, D and are asymptotically given by 



En 



En 



for n 



oo . 



(30) 
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TABLE I: Pressure dependence of measured quantities and Qo and calculated quantities #1,2,3 and l± at zero temperature 
and zero 3 He concentration. The ratio of the 3 He effective mass to its bare mass m%(P)/mz is taken from Ref. 25], which 
is obtained by fitting experimental data |26j| to the zero concentration. The BBP parameter at zero concentration Qq(P) is 
determined up to 10 atm by the experiment in Ref. |29|. ao(P) for 15, 20 atm (indicated with daggers) are read from Ref. |3C|. 
while they are at 6.0% concentration and not extrapolated to the zero concentration limit. The couplings (71,2,3 in the effective 
Lagrangian and l± in Ea. ll^ are calculated from the measured values with the use of Eqs. il'-MII and $A'A) . 



P [atm] 


mt(P)/m 3 


a (P) 


9l 


92 


93 


I- 


1+ 





2.18 


0.288 


0.541 


0.403 


0.809 


-2.34 


0.566 


5 


2.31 


0.242 


0.567 


0.260 


0.685 


-2.45 


0.485 


10 


2.44 


0.218 


0.590 


0.138 


0.581 


-2.57 


0.413 


15 


2.54 


0.172 1 " 


0.606 


0.0170 


0.472 


-2.65 


0.338 


20 


2.64 


0.148 1 " 


0.621 


-0.0820 


0.384 


-2.74 


0.278 



The criterion for an angular momentum I of 3 He atom 
in which bound 3 He states are available is simply given 
by 



= g 3 —w 2 - 2 gi ^-wl - I 2 > . 

7774 



(31) 



Using the definitions of the couplings (|20[1 and 77, — 

777,3 + 777774, the criterion can be rewritten as I- < I < 
with 



l± 



w 



m 3 



777,4 



y o/i4 7774 



(32) 



An infinite number of bound energy levels with their 
asymptotic ratio e ~ 27T / K appear for each integer I sat- 
isfying l_ < I < l + . We note that binding energies them- 
selves are not determined within the effective Lagrangian, 
because it can not access the vicinity of the vortex core 
where the magnitude of the condensate changes. 



C. Determination of parameters and results for 
3 He- 4 He mixture 

We should determine the free parameters of our theory 
7773 and dA/d[i4 from experiments. We adopt the 3 He 
effective mass 7773 at zero temperature from Re f. | 25j] . 
which is obtained by fitting experimental data |26( to 
the zero 3 He concentration. The ratio to the bare 3 He 
mass m^(P)/m3 is cited in Table [flat various pressures 
P = 0, 5, 10, 15, and 20 atm. 

On the other hand, the derivative of the 3 He energy 
cost with the He chemical potential dA/dp,4 is related to 
the relative fractional molar volume of 3 He in a 3 He- 4 He 
solution, or the Bardeen-Baym-Pines (BBP) parameter 
|27|. at zero concentration a by [2^ 



dA 

<9/J4 



a (P) + 1 . 



(33) 



ctQ has been determined as a function of the pressure P 
up to about 10 atm by the experiment pfij . and is shown 
in TableQJ The values for 15, 20 atm in TableHJare quoted 



TABLE II: Asymptotic ratios of energy levels Eoo/ -E00-1 = 
g-2-ir/K Q f DOunc i 3jj e s tates for each possible angular momen- 
tum at pressures P = 0, 5, 10, 15, and 20 atm. 



P [atm] 


I = 




I = -1 




/ = -2 







4.24 x 10" 


■3 


1.30 x 10" 


-2 


1.15 x 10" 


-3 


5 


3.13 x 10" 


-3 


1.38 x 10" 


-2 


2.62 x 10" 


-3 


10 


2.25 x 10" 


-3 


1.48 x 10" 


-2 


4.48 x 10" 


-3 


15 


1.32 x 10" 


-3 


1.46 x 10" 


-2 


6.11 x 10" 


-3 


20 


7.42 x 10" 


-4 


1.47 x 10" 


-2 


7.85 x 10" 


-3 



from Ref. [23|, while they are at 6.0% concentration and 
not extrapolated to the zero concentration limit. 

Then, l± in Eq. (|32|l as well as the couplings (71,2,3 in 
our effective Lagrangian (|20|l are determined from exper- 
imental values as functions of the pressure, which are 
listed in Table [I] We have put w = 1 because only 
\w\ = 1 vortex is energetically stable. The results show 
that bound 3 He states on a vortex of the 4 He supcrfluid 
are realized in / = 0, I = —1, and / = — 2 channels 
in a whole range of pressure, 0-20 atm, where experi- 
mental data are available. Note that because parity is 
broken by the vortex, the / = 1 and I = 2 states are not 
bound while the / = — 1 and I = —2 states are bound. 
The asymptotic ratios of energy levels (13011 of bound 3 He 
states for each possible angular momentum are shown at 
P = 0,5,10,15, and 20 atm in Table HU We conclude 
that those values are in principle measurable by accurate 
experiments. 



III. COLD FERMI GAS 

The effective field theory described above makes no 
assumption about the nature of the bosonic (tp) and 
fcrmionic ("0) degrees of freedom. This means we can use 
it for a two-component Fermi gas when tp is the phase 
of the Cooper pair and ip represents the fermion quasi- 
particles of a chosen component (say, the spin-up compo- 
nent if the two components correspond to different spins) . 
The fact that the Cooper pair is made up from the two 



fermions, one of which is the if) fermion, is of no impor- 
tance from the point of view of an effective field theory. 

Let us consider the cold two-component Fermi gas with 
a scattering length a, which has been recently archived by 
experiments using the technique of Fcshbach resonance 
d H E3, E1II1 Its ground state is found to be the 
superfluid in a whole range of a via usual BCS mechanism 
for a < (BCS regime) or Bose-Einstein condensation of 
tightly bound Coop er pairs (molecules) for a > (BEC 
regime) HHHHp. 

One can have in mind the situation of a slightly un- 
equal number densities in the two fermion components. 
Then, the ground state would be a homogeneous mix- 
ture of superfluid Cooper pairs and extra dilute fermions 
carrying the single component. Since the minimum of 
the fermion dispersion curve sits at zero momentum suf- 
ficiently deep in the BEC regime (corresponding to phase 
III in Rcf . 34] ) , such a system turns out to be within the 
scope of our low-energy effective field theory for boson- 
fcrmion mixtures. As in the case for 3 Hc- 4 Hc mixture, we 
are interested in the situation when the difference in the 
number densities is very small so that the interactions 
among the extra fermions are negligible. 

In our subsequent discussion, we refer to the inequality 
in the number densities as "polarization" according to the 
terminology of Ref. |34j in the case of spin-i fermions. 
Also, we shall use the term "BEC regime" for the case 
a > 0, and "BEC limit" for the limit na 3 -> +0. We 
introduce the notation \x = (/z-f + Mj.)/2 and H = (fi-\ — 
/i|)/2, with (itf and /xj being chemical potentials of each 
component of fermions for later use. 



BEC limit 



Microscopic description of the system 



It would be instructive to start with the microscopic 
description of the system in the BEC limit na 3 <C 1, 
where n is the fermion number density without polariza- 
tion. Since the system is dilute in this limit, the bound 
molecules can be regarded as pointlike bosons. There- 
fore, the dynamics of the molecules could be described 
by the following local Lagrangian [35|: 



4m 



Mb 



9b 



(* t *) 2 . (34) 



Here the field * represents the superfluid molecule which 
has its mass 2m with m being the fermion mass, /ib = 
2/.i + Eq is the chemical potential of molecule with Eq be- 
ing its binding energy. The coupling of its self-interaction 
gb is characterized in terms of the two-body scattering 
length between molecules by gb = 27rab/m. 

If we introduce the chemical potential for the polar- 
ization H larger than the energy gap of a single fermion 
\H\ > Eq/2, extra fermions carrying one sign of spin will 
be created on the top of the BEC ground state. The La- 
grangian describing such fermions and their interaction 




FIG. 1: (Color online) Solutions of Eq. 1381 as functions of 
r/£ for w = 1, 2, and 3. 



with the superfluid molecules will be given by 



2m 



3bf* t *V ,t V'- 

(35) 



%jj = f/;j or if>i depending on the sign of H . Hereafter 
H > is assumed to be positive without losing generality. 
The coupling of the interaction between the extra fermion 
and the bound molecule t/bf is characterized in terms of 
their two-body scattering length by gbf = 37rabf/m. Self- 
interactions among fermions are negligible in the dilute 
limit of extra fermions. We should note here that the 
mass of the extra fermion is provided by its bare mass 
m, because interaction effects become infinitely small in 
the BEC limit no 3 -> 0. 

The equation of motion from Eq. (|34J) leads to the 
Gross-Pitaevskii equation 



(i«9 +Mb)*(*,r) 



d 2 *(i,r) 
4m 



■3b|*(t,r)r*(t,r). 

(36) 



In order to consider a single vortex in the molecular su- 
perfluid, we set \E'(t,r) = y/nbh(r)e lw6 in cylindrical co- 
ordinates, rib = n/2 = fib /9b is the density of molecules 
far away from the vortex core and w is a winding number 
of the vortex. h(r) is some function of the radius r from 
the vortex core, which satisfies the boundary conditions 



h(r -► 0) -> 



and 



h(r 



1 



(37) 



Substituting *(t,r) into Eq. we obtain an equation 
for h[x), 



d 2 h I dh 

dx 2 x dx 



w 



h-h 3 



0. 



(38) 



where x = r/£ is the dimensionless radius normalized 
by the healing length of the molecular superfluid £ = 
l/y/4mg h n h . 
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TABLE III: Bound energy levels e n and their ratios e n /e n _i 
in the 5-wave channel (I = 0) for 7 = 1.47. 



n 




&ra/ £n-l 


u 


/i v in -1 
— 4.yo x iu 




1 


-3.39 x 10~ 3 


6.88 x 10~ 3 


2 


-1.91 x 1(T 5 


5.64 x 10~ 3 


3 


-1.08 x 10~ 7 


5.65 x 10~ 3 


4 


-6.10 x 10~ 10 


5.65 x 10~ 3 


5 


-3.45 x 10~ 12 


5.65 x 10" 3 


6 


-1.95 x 10" 14 


5.65 x 10~ 3 


00 




5.65 x 10" 3 



TABLE IV: Bound energy levels e„ and their ratios e n /e n _i 
in the P-wave channel (Z = ±1) for 7 = 1.47. 



n_ e n e„/e„_i 

"0 -1.57 x 10 -2 = 

1 -1.55 x 10~ 6 9.86 x 10~ 5 

2 -1.68 x 10~ 10 1.09 x 10~ 4 

3 -1.83 x 10~ 14 1.09 x 10~ 4 

00 1.09 x 10~ 4 



2. Bound fermion states on a superfluid vortex 



This equation should be solved under the boundary 
conditions JSJJ. Asymptotic forms of h(x) at x — > and 
x — > 00 are easily read from Eq. I)38|l as follows: 



h(x) oc x v 



for x^0, (39) 
for x — ► 00 . (40) 



The solutions of Eq. 138(1 connecting these two limits can 
be obtained numerically which are shown in Fig. ^ for 
w = 1,2, and 3. 

Because of the interaction between the molecule and 
the fermion, the vortex structure in the molecular super- 
fluid acts as a potential for the fermion. If the inter- 
action between the molecule and the fermion is repulsive 
(<7bf > 0), the fermion will be attracted to the vortex core 
because there are less molecules. Due to this effective at- 
traction, bound fermion states on the vortex turn out to 
be possible. The equation of motion from Eq. (|35l) with 
the solution of Eq. I|36|l gives the Schrodinger equation 
for the fermion under the vortex potential 



Eip(r) 



d 2 j>{r) 
2m 



+ g h{ n h {h(rf -l}^{r). (41) 



E is the energy eigenvalue of the fermion and is negative 
for bound states. 

Separation of variables and use of the normalized ra- 
dius x — r/£ lead to the equation for the radial wave 
function R(x) as follows: 



eR{x) 



dx 2 



ld_ 

x dx 



+ — +1 {h(x) 2 -l} 



R(x) 



(42) 



Here I is an angular momentum of the fermion and 
we have defined a dimensionless energy eigenvalue e = 
(2mE — k 2 ) £ 2 < and a coupling ratio 7 = gbi/2gh- k z 
is a momentum of the fermion along the vortex line. Note 
that Eq. (|42|l is invariant under This invariance 

cannot be exact because the vortex breaks parity. It is 
only an approximate property of the deep BEC limit, in 
which the fermion field feels only the magnitude of the 
condensate but not its phase [Eq. (|4ip]. 



The Schrodinger equation l|42|) with the asymptotic be- 
havior of the vortex potential l|40|) predicts an infinity 
number of energy levels for bound fermion states on a 
superfluid vortex at least in the S'-wave channel. We can 
rewrite Eq. (|42|l far away from the vortex core (x 3> 1) 
as 



d 2 R 

dx 2 



1BR 

x dx 



R = 



(43) 



with k = yj 7W 2 — I 2 . This equation has the same form 
as Eq. I|25|) so that the discussion in Sec. lIIBl is applicable 
here. Consequently, the criterion for the angular momen- 
tum I of the fermion in which bound fermion states are 
available is given by 



I 2 > 0. 



(44) 



An infinite number of energy levels for the bound fermion 
states appear with each integer I satisfying I 2 < jw 2 and 
their ratios are asymptotically given by 



= -27r/K 



for n — > 00 . 



(45) 



Whether we have bound fermion states with a certain 
I depends only on the value of the coupling ratio 7 = 
5bf/2gb = 3abf/4ab > 0. In particular, one can conclude 
that they are always possible for the S'-wave channel (I = 
0) in the BEC limit. 

Now we restrict ourselves to the energetically stable 
vortex k; = 1. Since the system is weakly coupled in 
the BEC limit na 3 — > 0, the scattering length between 
a molecule and a molecule or fermion is calculable as a 
function of the fermion scattering length a and results 
in Ob = 0.60-0.75 a [M |H H| or a b f = 1.179 a [H 
lifl l4l) . These values leads to the coupling ratio 7 = 
1.18-1.47, which predicts bound fermion states in the P- 
wave channel (I = ±1) as well as in the S'-wave channel. 
Hereafter we adopt the value a b = 0.60 a which has been 
confirmed by a quantum Monte Carlo simulation |42j| and 
experiments |9j. Then, the asymptotic ratios of bound 
energy levels eooAoo-i are given by 5.65 x 10~ 3 for the 
S-wave channel and 1.09 x 10~ 4 for the P -wave channel. 

The absolute values of binding energies should be de- 
termined by solving the Schrodinger equation (|42|l nu- 
merically. Results on bound energy levels for e„ < 10~ 15 
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FIG. 2: (Color online) Solutions of Eq. (M2> in the S'-wave channel (I = 0) for 7 = 1.47. Left panel: Bound energy levels e n 
for n — 0, 1, ... ,6 (dotted lines) and the potential energy 7 {h(xf - 1} (solid curve) are shown as functions of x = r/£. The 
normalized wave function of the ground state (n = 0) is also shown by the dashed curve. Right panel: \e„\ (dotted lines) are 
shown in the log scale. 




FIG. 3: (Color online) Solutions of Eq. g2J in the -P-wave channel [l — il) for 7 — 1.47. Left panel: Bound energy levels e n 
for n = 0, 1, 2, 3 (dotted lines) and the potential energy 1/x' 2 + 7 {h(x) 2 — 1} (solid curve) are shown as functions of x — r/£. 
The normalized wave function of the ground state (n = 0) is also shown by the dashed curve. Right panel: |e n | (dotted lines) 
are shown in the log scale. 



and their ratios e„/e„_i are listed in Table IIIII for the 
S-wave channel and in Table llVl for the P-wave channel. 
The ratios of bound energy levels e n /e n -i rapidly con- 
verge to their asymptotic values from n > 2 for both the 
S- and P-wave channels. Even in the first ratios ei/eo, 
their deviations from the asymptotic values are only 22% 
for the S'-wave channel and 9.1% for the P-wave channel. 

The bound energy levels (dotted lines) are shown in 
Fig. [3 for the S'-wave channel and in Fig. [3] for the P- 
wave channel as well as the potential energies I 2 jx 1 + 
■y {h(x) 2 — 1} (solid curves) as functions of x = r/£. 
While the binding energies of excited states (n > 1) 
seems degenerated into zero in the linear scale (left pan- 
els) for each channels, they are allocated with equal in- 
tervals 27t/(k1h10) in the log scale (right panels) as in- 
dicated by Eq. (|45|l . The wave functions of the ground 



state (n = 0) normalized as 

/>oo 

/ dxxR{xf = l (46) 

are also shown in each figure (dashed curves). In the 
S-wave channel, the ground state wave function is well 
localized around the vortex core r/£ < 5 and its binding 
energy is comparable to the bottom of potential energy 
I e o I / T = 0.334. Thus, we conclude that those bound 
fcrmion states in the BEC limit could be measurable by 
future experiments at least in the ground state of the 
S-wave channel. 

Away from the BEC limit, the fermion mass changes 
to the effective mass m*, which is larger than the bare 
mass to* > to [34|]. Accordingly, the coupling ratio in 
the Schrodinger equation l|42|l changes to the effective 
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one 7* = (m*/m)"/. Because 7* > 7, bound fermion 
states with higher angular momenta are expected to ap- 
pear as one departs from the BEC limit. However, if one 
is away from the BEC limit na 3 > 1, molecules overlap 
each other and then the microscopic description with lo- 
cal interactions Eqs. (|34(1 and (|35(1 will no longer be valid. 
In order to treat the problem away from the BEC limit, 
the effective field theory in terms of low-energy excita- 
tions is efficient. 



3. Connection to effective field theory 

Let us see how the effective field theory in terms 
of low-energy excitations, superfluid phonons and extra 
fermions, emerges from the microscopic Lagrangians <|34[) 
and l|35l) . For that purpose, we parametrize the field of 



bound molecule as ^(t,r) ~ \/nb(t, r)e 2ltptyt ' r ^ and as- 
sume that the magnitude 7ib(i, v) is slowly varying in 
time and space, which is satisfied far away from the vor- 
tex core. Then, Eqs. I(34|) and 1)35(1 can be rewritten by 



C-r + Ci 



(0<pY 



9b 2 



lift d$vb \&>p\* 



+(v + H)^ip. (47) 
I Im 

Since now the density of bound molecules rib is a varia- 
tional parameter, it should be determined by minimizing 
the action and results in 

(d(p\ 2 

g h n h = /i b - 2ip gbfp^ip ■ (48) 



in 

Substitution of this expression into Eq. I|47ll leads to 
an effective theory in terms of the superfluid phonon tp 
and the fermion excitation -0 as follows: 

H^-t - 

2 Y 2 



P( M ) + 4f ^ - J - [dtp? 



2m 



+ (p + H) ift$ (49) 



dA 



(dip) 2 
2m 



We have dropped a total derivative term and higher or- 
der terms in tp. Here P(fJb) = fj, 2 /2gt, is identical to the 
pressure of the molecular superfluid up to an irrelevant 
constant, and A(/x) = gbif-b/gb is the energy cost to in- 
troduce a single fermion into the superfluid originating in 
its interaction with the bound molecules. The low-energy 
parameters in the Lagrangian turn out to be given by 

o dn , n dP , . 

/ f 2 = — , 2 = - with n=— . 50 

o/j, m o/i 

This is the manifestation of general properties of Sec. Ill Al 
in the molecular superfluid case. The coupling between 
the superfluid phonons and extra fermions is also given by 
the derivative of the fermion energy cost A'(/i) = 2(?bf/5b 
as in Eq. (JT^J. 



B. Away from the BEC limit 

1. Effective field theory 

As mentioned above, the problem of writing down an 
effective field theory that couples the superfluid phonons 
(p and the extra fermions ip is identical to the same prob- 
lem in the 3 He- 4 He mixture. This is because their sym- 
metries and the pattern of symmetry breaking in the two 
cases are the same. Let us elaborate on this point in 
more detail. 

In the 3 He- 4 He mixed system, there are two conserved 
charges: the number of 4 He atoms N4 and the number of 
3 He atoms A3. The first charge is spontaneously broken 
by the 4 He superfluid ground state. The order parameter 
carries a unit N4 charge, but is neutral with respect to 
the A^3 charge. The 3 He atoms carry only the N3 charge 
which is unbroken, while being neutral with respect to 
the 4 He charge. 

In the slightly polarized fermion system, the number 
of spin-up fermions N-\ and the number of spin-down 
fermions N± are conserved separately. We suppose all 
extra fermions have spin up. It is, however, more conve- 
nient to use another basis, consisting of N± and the total 
polarization Y = A^ — N±. The Cooper pairs carry unit 
N± charge and its condensation spontaneously breaks this 
symmetry. The extra fermions carry a unit Y charge but 
are neutral with respect to N±. Thus, N± is equivalent to 
A^4 in the 3 He- 4 He mixture, and Y = Nf — N± is equiva- 
lent to A3. 

It is now possible to write down the effective La- 
grangian for the case at hand by direct analogy with the 
case for 3 He- 4 He mixture. We simply need to make the 
following replacement for the particle masses and chem- 
ical potentials: 



!/l4 



2m, 
2/i, 



m 3 



m, 

fL + H. 



(51) 
(52) 



The equations in Eq. (|52|l follow from the requirement 
that + fi3n3 — > /x-j-n-f + /ix n l when — > n\ and 

ri3 — > n-f — ri\. For convenience, we shall also make some 
rescalings so as to absorb extra factors 2 appearing in the 
equations: 
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2<P, 277 



(53) 



As a result, the Galilean invariant effective Lagrangian 
(fTHfl becomes 



r ft -2 



f 2 2 

+ (fi + H — A) iftip 
iip^dip 



iip^d ip \dip\ 



2m* 



(54) 



gx dip ■ 



2m 



.92 ip + .93 



(dipy 



2m 



iftip . 



Here m* is the effective fermion mass and the couplings 
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in the effective Lagrangian are given by 



with 



m d A d A 2 to , . 
5i = *7 — * ' 52 = -e V j 53 = t; »7 — , (55) 

with ?7 being defined by to* = (1 + rf)m. The universal 
relation among these couplings 



53 = 92 + 9i 



(56) 



is a consequence of the Galilean invariance. In particular, 
gi = and 52 = 93 = dA/dfi = 47 in the BEC limit 
m* — > to, which reproduce the result of Eq. Q49[l. 



2. Bound fermion states on a superfluid vortex 

Let us consider the bound fermion states on a super- 
fluid vortex with its winding number w. The discussion 
is parallel to that in Sec. Ill Bl What we have to be care- 
ful here is that the phase ip of the condensate around the 
vortex is given by tp(t,r) = w9/2 in cylindrical coordi- 
nates. Extra one-half compared to the 3 He- 4 He mixture 
case is because tp is normalized to be half of the phase of 
the Cooper pair iff = \^\e 2t(p in Eq. i|53|) for the fermion 
superfluid case. The equation of motion for the fermion 
field ip from Eq. (|5*4*Jl gives a Schrodinger equation as fol- 
lows: 



Etp(r) 



d 2 ip(r) igxw dip(r) g^w 2 



2m* 



2mr 2 d0 



8mr- 



■V(r). (57) 



The third term corresponds to the 1 jr 2 vortex potential 
far away from the vortex core in Eq. I|41l) in the BEC 
limit. On the other hand, the second term, which was 
absent in the BEC limit, is required by the Galilean in- 
variance when to* 7^ to. This term represents that the 
fermion with angular momentum in the opposite direc- 
tion to the vortex's one is energetically favored, because 
it decrease the total angular momentum of the system. 

Separation of variables leads to the equation for the 
radial wave function R(r) as follows: 



(2m* E - kfj R(r) 
d 2 



ld_ 

dr 2 r dr 



— 



5i- 



TO* wl 



53 



4to r 2 



(58) 
R(r). 



Here I is an angular momentum of the fermion and k z is 
its momentum along the vortex line. As we have proved 
in Sec. Ill Bl whether we have bound fermion states for a 
certain I is determined only by the sign of the coefficient 
of 1/r 2 . They are possible for angular momenta I which 
satisfy the criterion 



9 ^ W 



gi wl - I > 0. 

TO 



(59) 



Using the definitions of the couplings (|55|l and 77, to* = 
(1 + rj)m, the criterion can be rewritten as /_ < I < Z+ 



w 



1 to* 1 dAm* 

2 ~ 2m 2 V ~d[I~m 



(60) 



An infinite number of bound energy levels with their 
asymptotic ratio e - 27T / K appear for each integer I sat- 
isfying Z_ < I < l + . 



3. Discussion on parameters and conjectures 

We have two free parameters which cannot be de- 
termined within our effective field theory, the effective 
fermion mass to* /to and the derivative of fermion en- 
ergy cost A'(/x) in the superfluid. These two parameters 
are functions of the fermion scattering length a and the 
fermion chemical potential /1 or density n. They should 
be measurable by experiments in the same way as m| 
and «o = A'(/i4) — 1 in the 3 He- 4 He mixtures. While 
such measurements are not achieved yet, we still know 
some general properties of those quantities. 

As we have discussed, the effective fermion mass 
to*(/x, a) coincides with its bare mass to in the BEC limit 
na 3 <C 1. On the other hand, it has been argued that 
m*(jJ,,a) becomes infinite at the so-called splitting point 
(SP), located on the BEC half (a > 0) of the crossover 
diagram [34|. Therefore, it is natural to assume that 
to* (/1, a)/m is an increasing function of na 3 from unity 
(the BEC limit) to infinity (the SP limit), as na 3, in- 
creases from zero to some critical value. 

As for the derivative of fermion energy cost A'(/x, a), 
it is presumably positive in the BEC regime. This can 
be seen by writing the derivative of fermion energy cost 
as 



dA((i,a) dndA(n 7 a) 



d/J, 



dfi dn 



(61) 



Thermodynamic stability implies dn/d/i > 0. The en- 
ergy cost A to introduce a fermion into the superfluid 
originates in the interaction of the fermion with bound 
molecules in the BEC regime. Their interaction is consid- 
ered to be repulsive because the Pauli principle between 
the extra fermion and a fermion in the bound molecule 
with the same sign of spin acts as an effective repulsion. 
Thus, the fermion energy cost will be positive in the 
superfluid and an increasing function of the density n, 
which results in dA/dn > 0. In fact, this is the case 
in the BEC limit where A = gbfib with the repulsive 
coupling gbf = 1.179 a > 0. It is also true at the unitar- 
ity limit (na 3 = 00) where A ~ n 2 / 3 . Accordingly, we 
conclude that A' (fx) is positive in the BEC regime. 

With the use of those properties on the parameters, we 
can derive some interesting predictions from the formula 
(|60|l . Let us restrict ourselves to the energetically stable 
vortex w = 1. Because of l± — y/j — 1.21 in the BEC 
limit, the ,5- wave (I = 0) and P-wave (I — ±1) bound 
fermion states are possible. On the other hand, l± goes 
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FIG. 4: (Color online) 1+ and I- in Eq. (16011 as functions of 
m*/m for dA/dfi = 47. m*/m = 1 corresponds to the BEC 
limit, while m*/m — > 00 corresponds to the splitting point in 
Ref. |34|. The curves are shown for m*/m > 0. 

to negative infinity in the SP limit l± — * —00, because 
771* /m — * 00 and dA/dfi should be finite without phase 
transitions. Therefore, Z_ must get across integers less 
than —1 where new bound energy levels appear, and 
must get across integers not greater than 1 where exist- 
ing bound energy levels disappear between the BEC and 
SP limits. In particular, an infinite number of negative 
and large angular momenta become ready for the bound 
fermion states in the vicinity of the SP limit. 

For illustration, let us employ a somewhat unjustified 
ansatz, A'(/_t, a) = 47, which is an extrapolation from the 
weak coupling BEC limit. Using the value 7 = 1.47, l± 
in Eq. (|fc>U|) are evaluated as functions of to* /to in Fig.0] 
Bound fermion states are possible with angular momenta 
I which are integers between I- and ?+. The term pro- 
portional to g\ in Eq. I|54|) suppresses bound fermions 
with positive angular momenta, while it enhances bound 
fermions which have negative angular momenta in the 
opposite direction to the vortex's one. Fig. 21 shows that 
bound fermion states with positive angular momentum 
disappear away from the BEC limit, while those with 
arbitrary higher negative angular momenta become pos- 
sible as one approaches to the SP limit to*/to — > 00. 

The fermion effective masses where bound fermion 
states appear to*/to [/_=/ or disappear to* /to \i += i are es- 
timated for the first four angular momenta I in Table [V] 
Especially, the iS-wave bound fermion states which have 
the deepest binding energy in the BEC limit are found to 
disappear when (73 changes its sign at to*/to = 7.77. We 
emphasize that quantitative results in Fig.0|and Table IVl 
may not be reliable due to the ansatz we employed. 

Finally, it would be interesting to compare our results 
in the polarized Fermi gas with those in the unpolar- 
ized Fermi gas from the Bogoliubov-de Gennes approach 
[r7lll^ |. The results in Ref. |18|| show that bound fermion 
states on a vortex are possible only for I = 0,-1 at 
77a 3 = 1/37T 2 in the BEC regime, while bound fermion 
states become possible for all Z < in the unitarity limit 



TABLE V: 


Estimated 


values of the fermion 


effective mass 


where bound fermion states appear m*/m ;_ = 


and disappear 


m*/m ;+=i 


for several 


angular momenta I. 




I 


m* /in | ; 


1 


m*/m| ;+=i 


-2 


1.77 


1 


3.62 


-3 


2.88 





7.77 


-4 


4.09 


-1 


11.1 


-5 


5.37 


-2 


17.0 



77a 3 = 00. The tendency that possible angular momenta 
for bound states increase as one gets away from the BEC 
limit is consistent with what we found, however, we pre- 
dict I — 1 bound fermion states as well as I = 0, —1 in the 
BEC limit. How these two results match with each other 
as a function of the polarization will be an interesting 
future problem. 

IV. SUMMARY AND CONCLUSIONS 

In this paper, we constructed low-energy effective field 
theories for dilute fermion excitations in supcrfluids with 
emphasis on the Galilean invariance [Eq. (|18|l for boson 
superfluids or Eq. (|54"jl for molecular superfluids] . The 
Galilean invariance as well as the global symmetries of 
the system considerably restrict the possible form of the 
effective Lagrangian. We showed three terms represent- 
ing interactions between the fermion excitations tjj and 
superfluid phonons <p appear to the lowest nontrivial or- 
der of the fields. The couplings for the three interac- 
tion terms are not independent as a consequence of the 
Galilean invariance. They are written by the effective 
fermion mass to* /to and the derivative of the fermion 
energy cost with the chemical potential dA/d/j, [Eq. i|20|) 
or (|55|) ]. which are measurable quantities by experiments. 

We consider that the effective Lagrangian obtained 
here is valid as long as the following conditions are satis- 
fied, (i) The low-energy dynamics of superfluids is dom- 
inated by the excitations of superfluid phonons or, in 
other words, the variations in the magnitude of the con- 
densate are negligible, (ii) The fermion's dispersion is 
given by the quadratic power of its momentum with a 
positive coefficient as in Eq. (Q. Higher order correc- 
tions to the dispersion are negligible in the low-energy 
dynamics of fermions. (iii) Fermions are dilute or weakly 
coupled so that their self-interactions are negligible. 

With the use of the effective field theory, we studied 
bound fermion states on a single vortex of the superfluid. 
We derived a simple criterion for an angular momentum 
I of fermion in which bound fermion states are available 
[Eq. i|31|) or l|59|l ] . An infinite number of bound energy 
levels appear for angular momenta satisfying the criterion 
and their ratios are asymptotically given by Eq. I|30|l . 

We applied our effective field theory to two boson- 
fermion mixed systems, a dilute solution of 3 He in He 
superfluid and the BEC regime of a cold polarized Fermi 
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gas. For the 3 He- 4 He mixture, we determined param- 
eters of the effective theory from experimental data as 
functions of pressure. As a result, we predict that bound 
3 He states with I = —2,-1,0 will be realized on a vortex 
of the 4 He superfiuid in a whole range of pressure, 0-20 
atm, where experimental data are available. Asymptotic 
ratios of bound energy levels are calculated in Table |n] 
for each angular momentum. Those properties should be 
in principle confirmed by future accurate experiments. 

As for the cold polarized Fermi gas, we determined pa- 
rameters of the effective field theory from the microscopic 
description in the BEC limit. As a consequence, S-w&ve 
(I = 0) and P-wave (I — ±1) bound fermion states turned 
out to be realized in the BEC limit. Since the fermion 
mass compared to its bare mass m* /m could change from 
unity (the BEC limit) to infinity (the SP limit) in the 
BEC regime, bound fermion states with arbitrary higher 
negative angular momentum will become available away 
from the BEC limit. Especially, the bound fermion states 
with / = 0, ±1, which are realized in the BEC limit, are 
shown to disappear as one get away from the BEC limit. 

While we have concentrated on the study of bound 
fermion states on a superfiuid vortex in this paper, our 
effective field theory will be useful to study other prob- 
lems in the boson-fermion mixtures. Also, our effective 
field theory is universal to any other boson-fermion mixed 
systems. Those investigations should be performed in fu- 
ture works. 



rewrite Eq. IjA.ljl as 



a (P) 



lim 

x— >0 



V 3i {x,P)-V i {P) 
xV 4 (P) 



(A.2) 



and calculate V34 in the leading order of x. 

The chemical potential of the mixture is defined as 



H{x, P) = Xfl 3 (x, P) + (1 - x)^{x, P) 



(A.3) 



with /i3 and /i4 being the chemical potentials of the 3 He 
and 4 He atoms. We can calculate the molar volume of the 
mixture by differentiating the mixture chemical potential 
with the pressure as 



V 34 (x,P) = A 



dP 



(A.4) 



where A is the Avogadro's constant. With the use of the 
Gibbs-Duhcm relation 



dx 



(1 



dx 



(A.5) 



Eq. i|A.3|) is rewritten as 

ti{x, P) = (1 - z)M4(0, P) + xfi 3 {x, P) 

x' dfx s (x',P) 



(1 



1 ) / dx' - 
'0 1 



dx' 



(A.6) 
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APPENDIX: THE BARDEEN-BAYM-PINES 
PARAMETER 

Here let us remind of the relation between the energy 
cost A to introduce a 3 He quasiparticle into the 4 He su- 
perfiuid in Eq. l|12f> and the relative fractional molar vol- 
ume of a dilute solution of 3 He in 4 He, or the Bardeen- 
Baym-Pines (BBP) parameter, a |23,|23|. a is a directly 
measurable quantity by experiments. It is sufficient for 
the present purpose to limit ourselves at zero tempera- 
ture for simplicity. 

BBP parameter a(x, P) at the 3 He concentration x and 
the pressure P is defined through the following equation: 



V 34 {x,P) = V 4 (P) [l + xa(x,P)], 



(A.l) 



where V 34 and V 4 are the molar volume of the dilute 
solution of 3 He in 4 He and the pure 4 He, respectively. 
Total number of 3 He and 4 He atoms N is fixed here. 
Since we have interest in the dilute limit of a. let us 



E(k) = A(P) 



fc 2 



2m* 3 (P) 



■xE int (x,P,k) 



(A.7) 



The first term A is the energy cost to introduce a single 
3 He quasiparticle to the pure 4 He superfiuid, and the 
third term proportional to x represents the contribution 
of 3 He- 3 He interaction to the excitation energy. Since the 
third term gives higher order corrections of the order of 
0(x 2 ) to a(x, P), we neglect it hereafter. Therefore, the 
3 He chemical potential in the mixture is of the form: 



fi 3 (x,P) = A(P)+ Mf (z,P). 



(A.8) 



fiF is the chemical potential of a free Fermi gas of mass 
m 3 and density xN, which is simply given, at zero tem- 
perature, by 



(3tt 2 Y) 2 / 3 

Mx ' p) = ^4ipy x 



2/3 _ 



/i F (i,p). 



n2/3 



(A.9) 



Substituting Eqs. I|A.8|) and (|A.9|1 into Eq. I|A.6|) . we 

obtain 



H(x, P) = (1 - x)m(0, P) + xA(P) 

+ lx 5 '* fi F {l,P) + 0{x 2 ). 
5 



(A.10) 
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The differentiation of this equation with P leads to the 
expression of the molar volume of the mixture as follows: 

V 34 (x, P) = V 4 (P) - xV 4 (P) + xA^p- + 0(x 5 ^) , 



(A.ll) 



where V A (P) = Adfj, A (0, P)/dP. Consequently, BBP pa- 
rameter at zero 3 He concentration is given by 



a {P) 



A dA(P) 
Vi(P) OP 



- 1 = 



dA(P) 
d^(0,P) 



- 1 



(A.12) 



«o has been determined as a function of the pressure P 
by the experiment [29| up to about 10 atm, and found to 
be fitted very well by the following polynomial [Eq. (60) 
in [U]: 



a 



(p) = j2 a > pl 



(A.13) 



i=0 



with 



a i 
(is 

«4 



2.88069661 

-1.26990768 x 10 -2 atm" 1 



= 9.57212464 x 10~ 4 atm" 
= -6.13373794 x 10~ 5 atm 



(A.14) 



-3 



2.26173741 x 10" 6 atm" 4 



The unit is converted from "kgf/cm 2 " in 29] to "atm" 
here. The values of ao(P) for P = 0,5 and 10 atm are 
shown in Table [I] 



Finally, we note that A and /i4 in Eq. ljA.12fl are identi- 
cal to A and /14 introduced in Eq. (|12(l of the text. There- 
fore, the coupling between the 3 He quasiparticle and the 
4 He superfluid phonons dA/d/i is related with the BBP 
parameter at zero concentration ao by <9A = «o + 1 • 
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